Non-linear response of single-molecule magnets: 
field-tuned quantum-to-classical crossovers 
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Quantum nanomagnets can show a field dependence of the relaxation time very different from 
their classical counterparts, due to resonant tunneling via excited states (near the anisotropy barrier 
top). The relaxation time then shows minima at the resonant fields H n oc nD at which the levels 
at both sides of the barrier become degenerate (D is the anisotropy constant). We showed that 
in Mni2, near zero field, this yields a contribution to the nonlinear susceptibility that makes it 
qualitatively different from the classical curves [Phys. Rev. B 72, 224433 (2005)]. Here we extend 
the experimental study to finite dc fields showing how the bias can trigger the system to display those 
quantum nonlinear responses, near the resonant fields, while recovering an classical-like behaviour 
for fields between them. The analysis of the experiments is done with heuristic expressions derived 
from simple balance equations and calculations with a Pauli-type quantum master equation. 

PACS numbers: 75.50.Xx, 75.50.Tt, 75.45.+j, 75.40.Gb 



I. INTRODUCTION 

Single-molecule magnets are metal-organic clusters 
containing a magnetic core surrounded by a shell of or- 
ganic ligands which isolates the clusters from one another 
(see, e.g., Refs. 0). The most studied is Mni 2 , whose 
core contains eight Mn 3+ and four Mn 4+ ions strongly 
coupled via super-exchange interactions. This gives a 
ground state spin S — 10, while the large Jahn- Teller 
distortion on the Mn 3+ sites leads to a strong uniax- 
ial anisotropy. The energy levels have a bistable struc- 
ture e m ~ —D m (at zero field) with an energy barrier 
U — eo — es — 70 K to be overcome for the spin rever- 
sal. At low temperatures, these systems show the typical 
behaviors of superparamagnets, such as blocking or hys- 
teresis, yet at a much smaller scale of size. In addition, 
they form molecular crystals in which all molecules are 
nearly identical and, in the case of Mni2 acetate, have 
their anisotropy axes z parallel to the crystallographic c 
axis. 

These properties make molecular magnets model sys- 
tems to investigate whether quantum phenomena, like 
tunneling, survive in mesoscopic systems @. As is well- 
known, tunneling probabilities decrease exponentially 
with the height of the barrier to be tunneled through 
(height that grows with the system size) Q. At the 
same time, external perturbations can induce decoher- 
ence that degrades the quantum behavior 0, |f|. Fur- 
thermore, an external magnetic field H z detunes ener- 
getically the initial and final states for tunneling (i.e., 
those having +m and — (to + n) spin projections along 
z). Actually, many experiments have shown that tunnel- 
ing takes place at those fields where states of opposite 
orientation are degenerate, H n ~ nH\ (n = 0,1,2,... 
with Hi — 2g/iB-D — 4200 Oe in MnvA whereas it is 
suppressed for intermediate fields [E II • 

Such a resonant tunneling enables the spins to ap- 
proach faster their equilibrium state, giving rise to steps 



in the hysteresis loops around H z = H n |ij and to max- 
ima in the linear dynamical susceptibility xi- I 11 our 
previous work holhj, we found that resonant tunneling 
at H z = induces an extra contribution to the nonlinear 
response X3i making it larger (in magnitude) than the 
equilibrium one and having peaks reversed with respect 
to the classical predictions 0, ^| . In Refs. 0, [n| the 
dependence of the response on temperature, frequency, 
and orientation at zero field was studied. In this article 
we show how the tunneling contribution to the nonlin- 
ear response can be switched on and off by varying an 
external field, tuning and breaking successively the reso- 
nances. 



II. 



EXPERIMENTAL DETAILS 



Single crystals of Mni2 acetate were grown following 
the same procedure described in |Tlj . In order to increase 
the signal, the magnetic measurements were done on a 
collection of oriented and glued crystals. In our previous 
experiments [l(j,lll| we extracted the zero- field Xsi^) by 
fitting the dc field dependent x(w) to a parabola. Clearly, 
this method is not applicable to study how X3 depends 
on the external magnetic field itself. For this reason, in 
the present experiments we resorted to the more tradi- 
tional method of measuring non-linear susceptibilities by 
detecting the different harmonics X2(2w) and X3(3w) of 
the response. In the absence of bias field one had xi = 
(see below). A nonzero iJ z , however, makes xi the lead- 
ing nonlinear term, and we will mainly focus on it. 

The fields, both dc and ac, were applied parallel to 
the common anisotropy direction z of the clusters. We 
employed the susceptibility option of a commercial mul- 
tipurpose measuring platform (PPMS) which uses a con- 
ventional inductive method. It enables applying ac fields 
of amplitude ho < 17 Oe, and to selectively detect several 
harmonics of the exciting frequency lu/2tt < 10 kHz. To 
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FIG. 1: Illustration of the method employed for measuring 
the nonlinear susceptibilities (here T = 8 K, H z — 100 Oe, 
and ui/2n = 2kHz). X2(2w) and xz{3ui) are obtained, re- 
spectively, from the slope and the quadratic coefficient of the 
second m2{2uj)/ho and third rri3(3w) / 'ho harmonics of the out- 
put signal, measured as a function of the ac field amplitude 
h - 



separate the intrinsic nonlinear response of the sample 
from the possible contamination due to non-perfect har- 
monicity of the exciting ac coil, we measured the output 



signals m 2 {2uj)/h and m 3 (3uj)/ho at several ho- This 
gives the sought-for intrinsic contributions X2 and X3 as 
the terms proportional to ho and h\ respectively. An 
example of this procedure is shown in Fig. ^ 



III. 



RESULTS AND MODELIZATION 



The so-measured linear and nonlinear susceptibilities 
of Mni2 at T = 8 K are shown in Fig. [5] The frequency- 
dependent xi shows maxima near the resonant fields, 
where it approaches the equilibrium Xi q - We also dis- 
play (panel b) the second and third harmonic compo- 
nents measured at lu/2it = 2 kHz. (The high noise-to- 
signal ratio prevented from obtaining reliable \3 data for 
H z > 1 kOe.) The equilibrium nonlinear susceptibilities, 
also shown, were obtained differentiating Xi q , measured 
at the lowest frequency lu/2it = 1Hz. 

We clearly see that the magnitudes of the harmonics 
increase in the neighborhood of the resonant fields Ho = 
0, Hi and H 2 , where states of opposite S z are degenerate 
and the tunnel channels open. Besides, in contrast with 
the behavior of xii both \2 and X3 become, near Ho, 
larger than x^ and X;j q - Thus, when resonant tunnel 
sets on, the multi-harmonic response of these molecular 
clusters is enhanced. 



In order to understand these results we have derived simple expressions for the susceptibilities. This was done by 
solving, as in Ref. |13|. a system of balance equations for the net population of the two anisotropy potential wells: 

XiM = xTtT~- ' X2(2w) = xT-rr 4 x i q 7T~r- — 1 o- V M 

1 + iwr \ + 2ilut (1 + 1 lot) (1 + 2i lot ) 

X3(3^) = X° 3 q ^— -X°i q 



3ilut ^ i (l + iwr)(l + 3iwT) 



2iW cq 2i(^r') 2 



_ y e 1 V I (2) 

A2 (l + 2iwT)(l + 3icjr) Al (l+iwr)(l + 2iwr)(l + 3icjT) ' v ; 

The equilibrium X °k = {d k M z /dH*)/k\ are the derivatives of the magnetization curve, while r, t', and r" are the 
relaxation time and its corresponding field-derivatives (all evaluated at the working field H z ). At H z = we have 
Xa q = [since M Z (H Z ) = -M z (-H z )] as well as r 1 = [from t(H z ) = t{-H z )]. Then X2(2w) -> as H z -> 0, while 
in X3 the last two terms vanish. Therefore these equations extend the expressions of Ref. 131 to nonzero bias fields. 

I 



The first Eq. JJJ gives the ratio Xi/Xi* 5 < 1- Besides 
Xi depends, via the product lot, on how far the spins 
are from thermal equilibrium. By contrast, the nonlin- 
ear susceptibilities \2 and X3 include also terms depend- 
ing on t' and r", i.e., on how sensitive r is to changes 
of H z . As a result, the relaxation time does not sim- 
ply "renormalize frequency" , as occurs with xi > but it 
modifies the magnitudes of the nonlinear responses. This 
effect is missed in modelizations of the nonlinear suscep- 
tibility that fail to include the field derivatives of r |15| . 
As we discuss next, extending our arguments at H z = 
of |10L ITl| , it is this property that makes the quantum X2 



and X3 qualitatively different from the classical ones. 

According to Eq. 0), the relaxation time of the mag- 
netic clusters can be estimated from xi as Xi/^Xi where 
Xi and x'{ are the real and imaginary components of the 
first harmonic. This r is shown in the inset of Fig. [2] 
The data show minima at the resonant fields, in contrast 
with the monotonous behaviour of r in classical spins. At 
zero field and finite temperatures, Mni2 spins are able to 
tunnel between those excited magnetic states (m ~ 2-4) 
for which such process is not blocked by the internal bias 
caused by dipolar and hyperfine interactions jlfij . This 
results in a effective barrier reduced by a few magnetic 
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FIG. 2: Upper panel: linear susceptibility of Mni2 measured 
at T — 8 K versus a magnetic field applied parallel to the 
anisotropy axis, x, u/2n = 1Hz (^equilibrium); *, 500 Hz; 
A, 1 kHz; •, 2 kHz. Solid symbols, real parts; open symbols, 
imaginary parts. The inset shows the relaxation time r, as ob- 
tained from xi/<^XjL.[ see Eq. 0] as well as calculated for clas- 
sical spins (line) 14]. Lower panel, second harmonic suscep- 
tibility measured at the same temperature. • and o, ^j(2t^) 
and X2(2w) at 2 kHz; x, equilibrium X 7 = (dxT/dH z )/2. 
Inset: ., Xs(3w) at 2kHz; x, X T = (d 2 xT/^ H z )/6- The 
dotted vertical lines mark the resonant fields Hi ~ 4200 Oe 
and H 2 = 1H\. 



levels, say U ~ £±3 — £±s, so that the thermo-activated 
relaxation gets faster (SU ~ 4K). Tunneling is however 
suppressed as soon as the external bias £ TO = 2g^mH z 
exceeds the tunnel splitting A m , slowing down the relax- 
ation (the full barrier has to be overcome). As a result, r 
is minimum at zero field, whence r" > 0, while t' changes 
sign from < to > 0. 

The same features are repeated every time the field 
brings magnetic levels again into resonance H n = n x 
(2<7/ib)-D. Therefore, tunneling becomes, at any crossing 
field, an additional source of nonlinear response via r 1 
and t". Besides, accounting for the signs of the r deriva- 
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FIG. 3: Theoretical calculations of X2(2co) [Eq. Jy] for quan- 
tum and classical spins. In the latter case, we used Brown's 
classical formula for r [l4( and in the former r was calculated 
by solving Pauli's quantum master equation, as described in 
|16j . Notice the difference in the %2 axis scale (the solid equi- 
librium curve is the same in both panels). 



fives and Eqs. Q and one sees that the sign of the 
nonlinear susceptibilities can be reversed with respect to 
the classical ones. (In the classical model r decreases 
monotonically with increasing field [l4j|. inset of Fig. |3 
giving t" < and r' < for any H z ; the same occurs 
in a quantum thermoactivation model not including the 
possibility of tunneling p^lisj]). 

It is interesting that both behaviors can be obtained 
in our case just varying the external field. For fields 
between resonances, tunneling becomes blocked for all 
states and the spins reverse by thermal activation over 
the total ( "classical" ) energy barrier. But when a cross- 
ing field is approached, the strong nonlinearity of r shows 
up with its characteristic contribution to the nonlinear 
susceptibilities via t' and r". 

To confirm this interpretation we have computed the 
nonlinear responses from Eqs. Q and © but incorpo- 
rating the relaxation time obtained by solving a Pauli 
quantum master equation (as in Ref. Il6|) . The results 
(Fig. 0) show that the quantum contribution to ^2(2^) 
is dominant near the resonances. This is due to the small- 
ness of the tunnel splitting of the relevant states for our 
Mn i2 sample: A 4 ~2x 1CT 2 K and A 2 ~7x 10 _1 K. 
This means that the fields required to block tunnel via 
these levels, albeit relatively small (~ 20 Oe and 1000 Oe, 
respectively), give rise to relatively large changes in r, 
and hence, large r' and t". In the "classical" regime, 



4 



by contrast, the scale of change of r is determined by 
the anisotropy field H a ~ (25 — As this is very 

large in Mni2 (~ 95kOe), one has a comparatively slow 
decrease of r with H z . This correspond to small values 
of the derivatives r' and t" and in turn of the classical 
(non-tunnel) nonlinear susceptibilities. 

It is also worth mentioning that the tunnel splittings, 
which determine the width of the r vs H z resonances, 
are further broadened by dipolar and hyperfine interac- 
tions 0| . In fact, the master-equation calculations tell 
us that tunneling via lower lying states would give rise 
to enormous spikes in \2 (Aio ~ 7 x 10 _10 K for the 
ground state ). However, these peaks are also very frag- 
ile, being easily suppressed by environmental bias fields 
and therefore not observed experimentally. 

IV. SUMMARY AND CONCLUSIONS 

We have studied experimentally the nonlinear suscep- 
tibilities of a Mni2 acetate molecular magnet in the pres- 
ence of a longitudinal field H z . The standard method of 
measuring the harmonics of the response to an oscillating 
field ho cos(uit) has been employed. By using several am- 
plitudes ho we managed to isolate the genuine nonlinear 
susceptibilities \2 and X3 (oscillating with 2ut and Scot) . 
The low signal-to-noise ratio (in spite of gluing several 
single crystals) prevented from obtaining good X3 data; 
fortunately, we obtained nice curves for \2, which is the 



leading nonlinear term when a bias field is applied. 

The analysis and interpretation of the susceptibility 
curves was done with help from expressions derived with 
a simple system of balance equations (for the potential 
well populations). At variance with previous formulae, 
the field derivatives of the magnetic relaxation time are 
captured by our expressions. This, together with the 
known strong effect on r(H z ) of resonant tunneling near 
the barrier top, permitted to understand the experimen- 
tal phenomenology. We also plugged in those equations 
the t(H z ) obtained solving a Pauli quantum master equa- 
tion for Mni2, supporting this interpretation. 

Near the resonant fields H n (matching the levels at 
both wells) the \2 vs. H z curves neatly amplify the res- 
onant tunnel, as this entails large dr/dH z . For fields in 
between the H n 's, tunnel is blocked and the response is 
governed by the thermo-activation over the total barrier, 
as in the classical case. This does not give such large 
t', while its sign is reversed with respect to the tunnel 
contribution. Thus the sensitivity of \2 to the local fea- 
tures of the t(H z ) curve provides an alternative method 
to asses if tunneling plays a role in the relaxation of a 
superparamagnet; and if so, in which field ranges it does 
take place. 
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